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SUMMARY 


Developing  a  highly  accurate  method  of  predicting  the  arrival 
time  of  a  LORAN  pulse  traveling  as  a  ground  wave,  this  report  considers 
the  one-dimensional  integral-equation  method  of  analyzing  propagation 
over  irregular  terrain.  Numerical  calculations  and  theoretical  devel¬ 
opments  demonstrate  how  topography  and  various  approximations  affect 
the  calculated  propagation  of  a  ground  wave.  We  find  that  curvature 
effects  entering  through  geometric  terms  of  the  integral  equation  for 
the  propagation  of  phase  work  strongly,  but  that  those  entering  through 
a  modified  impedance  condition  are  insignificant.  Dispersion  of  the 
pulse  can  cause  cycle  ambiguity  in  the  signal.  For  a  ground  conduc¬ 
tivity  of  0.001  mho/m,  cycle  ambiguity  appears  at  a  distance  from  the 
transmitter  of  roughly  1500  km. 

Extensive  numerical  results  are  presented  for  generic  terrain 
types  selected  to  illustrate  salient  features  of  solutions  to  the 
integral  wave  equation.  Results  for  simple  hills  show  that  arc-length 
correction  of  the  propagation  path  is  insignificant  compared  with 
geometric  effects  on  phase  speed  in  predicting  phase  retardation. 
Isolated  terrain  features  cause  large  local  perturbations  in  phase 
arrival;  far  beyond  the  irregularity  on  the  propagation  path,  effects 
are  observable  but  diminish  with  distance.  The  one- dimensional  ver¬ 
sion  of  the  integral  wave  equation  can  overstate  the  effect  of  terrain 
features  that  are  localized  in  the  direction  transverse  to  the  great- 
circle  propagation  path. 
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PREFACE 


This  report  continues  Pacific-Sierra  Research  Corporations 
(PSR)  analysis  of  ground  wave  propagation,  specifically  the  effort 
to  develop  a  highly  accurate  method  of  predicting  LORAN  coordinates 
over  irregular  terrain.  In  an  earlier  report.  Field  and  Allen  de¬ 
termined  which  approximations  in  the  chosen  mathematical  treatment 
have  an  important  effect  on  accuracy,  and  discussed  their  significance. 

The  present  report  continues  the  analysis  of  Field  and  Allen, 
and  determines  the  numerical  consequences  of  the  approximations  they 
treated.  In  three  essentially  independent  analyses,  this  report 
evaluates  the  importance  of  the  effect  of  ground  curvature  on  ground 
impedance;  defines  the  distances  over  which  group  effects  modify  the 
arrival  time  of  a  LORAN  pulse;  and  presents  numerical  results  ob¬ 
tained  from  PSR's  developing  propagation  code,  which  will  handle 
arbitrary  ground  data  and  has  been  completed  for  the  case  of  termi¬ 
nals  on  the  ground.  The  PSR  ground-wave-propagation  code  can  be 
incorporated  into  a  LORAN  map-making  program.  An  analysis  of  the 
continuity  of  solutions  for  ground-based  and  low-elevation  receivers 
is  included  as  an  appendix;  numerical  methods  for  dealing  with  ele¬ 
vated  receivers  will  be  addressed  in  a  later  phase  of  this  analysis. 


E.  C.  Field  and  R.  Allen,  Propagation  of  the  Low-Frequency 
Groundwave  over  Nonuniform  Terrain ,  Rome  Air  Development  Center, 
RADC-TR-78-68 ,  March  1978. 
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I.  INTRODUCTION 


LORAN  navigation  depends  on  the  determination  of  arrival  times  of 
pulses  centered  in  frequency  at  100  kHz.  Electrical  and  topographical 
variations  along  the  propagation  path  can,  however,  considerably  in¬ 
fluence  such  arrival  times.  Pacific-Sierra  Research  Corporation  (PSR) 
is  developing  a  numerically  efficient  method  of  predicting  ground  wave 
propagation  over  irregular  terrain  that  most  nearly  satisfies  two 
conflicting  goals:  (1)  achieving  the  highest  precision  commensurate 
with  the  2.5  nsec  resolution  of  modern  receivers  (Fehlner  et  al. 

[1973])  while  (2)  incorporating  general  electrical  and  topographical 
ground  characteristics. 

Virtually  exact  analyses  reported  previously  (e.g.,  Norton  [1937]) 
are  valid  only  for  flat  or  smooth,  spherical  surfaces.  The  one¬ 
dimensional  integral-equation  approach  used  here  accommodates  a  great 
variety  of  conductivity  and  topographical  features,  provided  they  are 
not  too  abrupt.  One  version  of  this  method  lias  been  explored  by 
Johler  and  Berry  [1967],  Our  approach  uses  different  numerical  tech¬ 
niques  and  considers  in  greater  depth  several  important  aspects  of 
LORAN  pulse  propagation  for  ground-based  terminals.  The  integral 
equation  approach  involves  several  approximations;  however,  careful 
handling  of  certain  terras  permits  reasonable  accuracy,  as  discussed 
by  Field  and  Allen  [1978]. 

This  report  extends  the  analysis  of  Field  and  Allen,  and  deter¬ 
mines  the  numerical  consequences  of  the  approximations  they  treated. 
Section  II  examines  the  impedance  boundary  condition.  Previously, 
the  only  definitive  treatment  of  curvature  effects  on  ground  impedance 
was  that  of  Rytov  [1940],  Unfortunately,  his  analysis  is  available 
only  in  Russian;  furthermore,  application  of  Rytov fs  formulas  requires 
ground  surface  data  to  be  in  an  impractical  form — given  by  the  metric 
in  an  arbitrary  set  of  curvilinear  coordinates  molded  to  the  surface. 
Thus  Rytov  obfuscates  physical  interpretation  of  his  conclusion.  We 
have  reexamined  the  problem  using  the  same  physical  assumptions,  sub¬ 
stituting  an  equivalent  but  more  practical  description  of  the  surface. 


cc^Lion  III  examines  dispersion  of  the  multiple-component  LORAN 
pulse,  which  cannot  be  treated  exactly  by  the  ground-wave  integral 
equation  useful  for  monochromatic  waves.  In  the  absence  of  disper¬ 
sion,  that  equation  can  be  applied  directly  to  the  propagation  of  the 
LORAN  pulse.  In  the  presence  of  dispersion,  however,  it  can  only 
approximately  predict  the  time  of  arrival  of  some  physical  property 
of  the  pulse — e.g.,  a  zero  crossing — because  such  properties  can  no 
longer  be  exactly  identified  with  phase.  If  the  pulse  were  nearly 
monochromatic,  the  approximation  would  be  very  good. 

LORAN  receivers  typically  use  a  particular  zero  crossing  as  the 
arrival  time  of  the  pulse.  However,  calculations  of  arrival  time  are 
usually  based  on  carrier  phase.  Our  analysis  shows  that  a  LORAN  pulse 
is  sufficiently  monochromatic  that  a  specified  zero  crossing  can  be 
closely  approximated  by  a  specified  point  of  the  carrier  phase,  pro¬ 
vided  the  range  is  not  too  great.  Nevertheless,  since  we  find  that 
zero  crossings  traverse  the  pulse  envelope,  there  is  a  range  beyond 
which  the  number  of  a  particular  zero  crossing  is  ambiguous. 

Section  IV  presents  numerical  results  for  ground-based  terminals. 
To  illustrate  the  effect  of  generic  terrain  features  on  the  phase  of 
a  100  kHz  sinusoidal  signal,  we  consider  a  succession  of  models  be¬ 
ginning  with  smooth  ground,  then  proceeding  to  increasingly  compli¬ 
cated  structures — i.e.,  a  flat,  sloping  ramp  inclined  with  respect  to 
smooth  ground;  a  ramp  whose  sections  have  different  slopes;  and  a 
pyramid-shaped  hill  consisting  of  two  equal-length  ramps.  Simple 
sets  of  ramps  can  be  constructed  to  yield  virtually  the  same  phase 
behavior  caused  by  any  more  complicated  ground  shape.  For  example, 
the  pyramid’s  effect  on  phase  has  a  signature  typical  of  any  simple 
hill  of  similar  height  and  width. 

Section  IV  also  considers  combinations  of  hills  and  valleys  on 
a  round  earth,  varying  their  size,  shape,  number,  spacing,  and  loca¬ 
tion.  The  numerical  importance  of  certain  mathematical  terms  is 
evaluated . 

Appendixes  A,  B,  and  C  present  derivations  of  the  curvature  im¬ 
pedance,  the  LORAN  spectral  distribution,  and  the  analytical  continuity 
of  he i glit  gain  at  the  ground  surface,  respectively. 
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II.  CURVATURE  CORRECTIONS  TO  SURFACE  IMPEDANCE  CONDITION 


In  discussing  curvature  terms  in  the  surface  impedance  formula, 
previous  wave  propagation  analyses  [Senior,  1961;  Leontovich,  1944] 
refer  ultimately  to  the  calculation  of  Rytov  [1940].  Rytov*s  calcu¬ 
lation,  though  quite  elegant,  is  not  easily  applied  to  real  terrain 
since  it  assumes  that  the  metric  is  given  for  the  surface  over  which 
propagation  is  considered.  Curvilinear  coordinates  are  then  employed 
to  calculate  a  curvature  correction  from  derivatives  of  the  logarithm 
of  metric  component  ratios. 

Rytov1 s  method  simplifies  the  problem  in  the  abstract,  but  com¬ 
plicates  the  solution  of  concrete  cases.  However,  using  a  straight¬ 
forward  formalism,  we  can  find  a  more  practical  solution  to  the  same 
problem.  To  do  so,  we  rederive  curvature  corrections  using  cartesian 
coordinates,  MRS  units,  and  the  following  physical  assumptions  em¬ 
ployed  by  Rytov: 

•  We  can  neglect  the  displacement  current  in  the  medium  over 
which  the  waves  are  propagating. 

•  The  index  of  refraction  is  large  enough  that  Snell* s  law 
gives  propagation  nearly  normal  to  the  surface  of  the  con¬ 
ductor. 

•  The  skin  depth  in  the  medium  is  small  relative  to  character¬ 
istic  radii  of  curvature,  so  we  may  expand  the  electric 
field  in  a  power  series  of  the  ratio  of  skin  depth  to  radii 
of  curvature. 

We  outline  here  the  method  of  calculating  the  first  curvature 
correction  to  the  surface  impedance  condition;  Appendix  A  contains 
a  detailed  derivation.  The  calculation  begins  by  solving  the  follow¬ 
ing  two  Maxwell  electromagnetic  field  equations: 
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V  x  E  =  -y 


3h 

dt 


(1) 


V  X  H  =  OE  , 


(2) 


where,  expanding  the  electric  field  E  and  the  magnetic  intensity  H, 


E  =  (Eq  +  SEj  +  62E2  +  •••)  exp  j  i  |  >  (3) 

H  =  (Hq  +  6Hj  +  62H2  +  •••)  exp  j  i  —  -  cut  J  j  ,  (4) 

and  the  skin  depth  6  is  given  by 


and  where  O  is  ground  conductivity,  measured  in  mhos  per  meter;  p 

is  permeability;  t  is  time;  ^  is  the  phase;  and  (a)  is  the  frequency. 

The  condition  of  normal  propagation  is  applied  by  defining  three 
~y  -+  ~~y  — 

vectors.  A,  B,  and  n,  where  A  and  B  are  orthonormal  vectors  tangent 

— V 

to  the  conductor  and  n  is  the  unit  normal  into  the  conductor  so  that 

— V  — ► 

A  x  B  =  -n.  Then  we  require  that  V\j)  x  n  =  0.  Figure  1  illustrates 
the  geometry  of  the  vectors. 

The  zeroth-order  terms  in  the  expansion  yield  -  0.  The  first 
order  terms  give  the  following  conditions: 


*1  •  0  -  *0  •  S  •  eia- -i/I7Jhob  -  eib'i,/I7Ihoz-  <6a> 


and 


V\p  •  n  *  ^2i  .  (6b) 

The  curvature  terms  enter  at  the  second  order  as  tangential  com¬ 
ponents  of  V  x  Hq.  We  may  thus  evaluate  the  electric  field  and 
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x 


Fig.  1 --Geometry  for  calculation  of  surface  impedance 
with  curvature  correction 


approximate  it  by  its  first  three  terms  (the  zeroth  term  vanishes); 
the  resulting  equation  is  the  analog  of  Rytov’s  Eq .  (12a): 

Ea  =  e1’*'  wyH0B  j-  i/C/2  6  +  <S2  £-i/T72  ^  |  (A  •  V  x  B  +  B  •  V  •  A) 

(7) 

Setting  the  curvature  of  the  incident  wave  to  zero  and  evaluating 
the  vector  products  in  Eq.  (7)  for  a  one-dimensional  hill,  we  arrive 
at 

£  •  ^[,+!  ■■•<»>]  ■  <8> 

where  y^  and  Cq  are,  respectively,  the  permeability  and  permittivity 

of  free  space;  /yn/en  is  the  impedance  of  free  space;  n  is  the  index 
u  U  g 

of  refraction  at  the  ground,  which  with  the  present  assumptions  is 


and  2M(y)  is  an  approximat ion  to  the  radius  of  curvature  for 

a  hill  with  gentle  slope,  where  the  elevation  function  of  the  hill  is 
* 

given  as  z(y). 

Equation  (7)  may  be  interpreted  in  terms  of  the  focusing  of  the 

wave  within  the  conducting  medium-  If  the  incident  wave  is  vertically 

polarized,  as  occurs  in  LORAN  transmission,  it  is  appropriate  to  use 

the  version  of  Eq.  (8)  that  relates  E  to  HA  ,  since  the  sign  of  z" 

d  UA  ^ 

appearing  in  Eq-  (8)  uses  the  convention  that  A  is  the  direction  of 
the  axis  of  the  hill.  The  appropriate  relation  is  found  by  using  the 
transformation  A  B,  B  +  -A,  E  -+  H,  H  -E.  Then 


(9) 


E_ 

B 

H 

OA 

f  [l  -  T  *"<»>] 

g  L  J 


This  expression  gives  an  easily  programmable,  modified  impedance  equal 
to  [1  -  (6zn/4)]  times  the  flat  surface  value.  Our  correction  factor 
is  slightly  more  accurate  than  that  given  by  Leontovich  [1944]  without 
derivation. 

Equation  (9)  may  be  interpreted  physically  in  terms  of  the  stretch¬ 
ing  of  the  electric  field  lines  at  the  surface,  shown  schematically  in 
Fig.  2.  The  field  is  unaffected  by  the  surface  curvature.  The 
field  in  the  ground,  however,  gets  stretched  (squeezed)  as  the  wave 
front  begins  to  diverge  (converge).  Energy  conservation  requires  the 
electric  field  to  change.  The  changes  in  the  ratio  E/H  may  thus  be 
anticipated  depending  on  whether  the  hill  acts  as  a  converging  lens 
(zM  <  1)  or  a  diverging  lens  (zM  >  1).  The  limit  of  validity  for  the 
series  expansion  is  that  the  wave  must  be  severely  attenuated  before 
ravs  <ross  in  the  ground  and  nonlinearities  enter  the  problem. 

The  discussion  in  Sec.  IV  uses  the  convention  z  as  a  function 
ci  x,  the  distance  on  a  line  from  the  transmitter  to  the  receiver. 
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Fig.  2--Physical  interpretation  of  surface  impedance  correcti 
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ll  I.  TIME  OF  ARRIVAL  OF  THIRD-CYCLE  ZERO  CROSSING 


For  accurate  LORAN  timekeeping  and  navigation,  the  time  of 
arrival  of  the  third-cycle  zero  crossing  of  a  LORAN  pulse  must  be 
predicted.  The  Hufford  equation  [Hufford,  1932],  used  here  (Sec.  IV) 
and  by  others  (e.g.,  Johler  and  Berry  [1967]),  yields  the  phase  cor¬ 
rection  for  a  purely  sinusoidal  signal.  A  LORAN  pulse,  however, 
necessarily  is  a  superposition  of  many  frequencies,  so  it  is  appro¬ 
priate  to  examine  how  accurately  sole  reliance  on  the  Hufford  equation 
at  the  LORAN  carrier  frequency  predicts  the  propagation  of  a  zero 
cross ing. 

* 

A  LORAN  pulse  F ( t )  can  be  approximated  by 
F(t)  *  sin  (u>t  +  0  +  0)  I  -  e 

L  p 

3 

where  oj  *  2tt  *  10  kHz, 

0  *  absolute  envelope  to  cvcle  difference, 

9  *  0  or  ii, 

*  time  to  peak  of  pulse  psec. 

Since  such  a  pulse  contains  significant  energy,  distributed  from  90 
to  110  kHz,  dispersion  of  the  pulse  must  be  considered  in  order  to 
accurately  predict  the  third-cycle  zero  crossing. 

A  propagating  Oaussian  pulse  may  be  defined  by 


o-'Vl 


where 
v  1  s 


0(xf  t )  *  exp 

k  anti  ii>  *re  t 
t  he  velocity  ol 


he  central  wave  number  and  frequency, 
the  pulse;  and  Ax  is  the  width  of  the 


respectively; 
pulse.  An 


Measurements  bv  the  Applied  Physics  Laboratory  (APL) 
et  al.,  19  76)  indicate  that  u>  *  u>(t)  and  is  not  constant, 
quency  drift  is  probably  due  to  the  particular  transmitter 
API/ s  experiment,  so  Is  not  included  in  our  analysis. 


[ McCartv 
That  f re¬ 
used  in 
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analysis  of  such  a  pulse,  familiar  from  discussions  of  quantum  mech¬ 
anics  (e.g.,  Bohm  [1951]),  shows  that  the  width  of  the  pulse  spreads 
according  to 


The  shape  of  the  LORAN  pulse  does  not  significantly  differ  from  that 
of  a  Gaussian,  so  our  formula  should  fairly  accurately  describe  the 
spreading  of  the  pulse  as  it  propagates. 

If  dissipation  is  ignored  or  is  independent  of  wave  number,  then 
the  center  frequency  is  unchanged  and  the  envelope  merely  spreads 
rather  than  shifting  center — thus  including  more  cycles  of  the  carrier. 
Furthermore,  the  center  of  the  pulse  moves  with  the  group  velocity, 
whereas  the  center  frequency  component  moves  with  the  phase  velocity. 
Since  the  two  velocities  are  not  necessarily  equal,  what  emerges  is  a 
spreading  envelope  moving  with  a  varying  time  (or  distance)  from  the 
front  of  the  pulse  to  any  chosen  zero  crossing. 

The  time  of  arrival  of  the  third-cycle  zero  crossing  will  be  the 
sum  of  the  time  of  arrival  of  the  front  of  the  pulse  (when  the  signal 
exceeds  the  threshold  of  detection)  plus  the  time  for  three  cycles  of 
the  carrier  frequency  to  pass.  That  sum  does  not  necessarily  deter¬ 
mine  the  time  of  crossing  of  the  third-cycle  zero  generated  by  the 
transmitting  antenna.  Since  cycles  are  being  added  to  the  pulse,  and 
the  zeros  traverse  the  pulse,  the  received  third-cycle  zero  need  not 
follow  the  front  of  the  pulse  by  the  same  time  interval  as  was  origin¬ 
ally  generated — indeed,  the  received  third-cycle  zero  could  be  different 
from  the  original  tagged  point.  Figure  3  illustrates  the  relationship 
between  the  times  of  arrival  of  the  third-cycle  zero  crossing  t^  and 
the  front  of  the  pulse  t^.. 


METHOD  OF  ANALYSIS 

To  analyze  the  above  zero-crossing  ambiguities,  we  must  know  the 
phase  and  group  velocities  as  well  as  the  spreading  factor  for  the 
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pulse.  All  are,  in  general,  obtainable  from  the  dispersion  relation 
for  a  wave.  Unfortunately,  the  Hufford  equation  is  not  a  dispersion 
relation  since  it  does  not  give  ca  as  an  algebraic  function  of  k. 
Therefore,  we  must  numerically  generate  a  dispersion  relation  by  com¬ 
puting  the  phase  of  the  signal  from  the  Hufford  equation. 

The  wave  number  k  of  a  wave  is  defined  locally  as  the  gradient 
of  the  phase 


where  the  field  $  is  equal  to  |$|e^.  The  Hufford  equation  yields  W, 
which  is  defined  by  $  *  2W$q,  and  may  considered  the  attenuation 
function  of  the  field  The  free  space  solution  is  Thus  if 

W  =  |w|e^,  then 


k(x)  =  M  +  w 

W  3x  c  • 

where  c  is  the  speed  of  light  in  a  vacuum. 

We  may  evaluate  the  function  k(x)  at  90,  95,  100,  105,  and  110 
kHz  and  fit  it  to  a  second-degree  polynomial  at  several  different 
distances  to  obtain  local  dispersion  relations  of  the  form 

a)  =  4-  C^k  +  C^k  , 

where  C^,  C2*  and  are  polynomial  coefficients.  The  group  velocity 

V  is  then 

g 


3u) 

3k 


100  kHz 


2C3k  +  C2  • 


the  phase  velocity  V 


P 


is 


100  kHz 


ci 

T  +  C2  +  c3k 
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2  2 

and  8  U)/3k  ■  2C^,  which  is  needed  to  determine  the  spreading  factor. 

Tables  1  and  2  list  the  resulting  quantities  for  two  different  ground 

conductivities  at  various  distances:  the  spreading  factor  in  units 

of  Ax/Axrt  -  1  and  the  computed  values  of  V  -  V  . 

0  K  g  P 

The  quantities  in  Table  1  show  that  pulse  spreading  is  not  signif¬ 
icant  over  the  propagation  distances  considered.  This  conclusion 
differs  from  that  of  Johler  and  Walters  [1959],  whose  calculation  showed 
more  pulse  spreading.  The  difference  arises  because  we  consider  here 
an  actual  LORAN  pulse,  whereas  they  considered  an  exponentially  damped 
cosine  signal  of  only  about  one  cycle  in  duration.  Larger  spread  is 
prevented  by  the  greater  length  of  a  LORAN  pulse,  as  shown  by  the 
(Axq)^  term  appearing  in  Eq .  (10).  It  is  then  accurate  to  picture  the 
pulse  as  an  envelope  propagating  along  with  zero  crossings  of  the 
carrier  signal  receding  (since  >  V^)  from  the  front  of  the  pulse. 

The  speed  with  which  the  phase  recedes  varies  with  distance 
propagated.  We  can  estimate  the  distance  at  which  phase  identity 
difficulties  are  introduced  into  the  problem  by 


For  a  ground  conductivity  of  0.001  mho/m,  this  equation  yields  a  dis¬ 
tance  of  1500  to  2000  km.  For  0  -  0.01  mho/m,  however,  numerical 
results  show  the  distance  to  be  about  ten  times  as  far. 

This  argument  is  strongly  suggestive  but  not  numerically  precise 
for  the  lower  conductivity  case,  because  such  concepts  as  phase  and 
group  velocity  are  approximate  in  nature  except  when  applied  to  pulses 
of  infinitely  narrow  spectral  width.  The  dispersion  relation  we  used 
was  found  numerically  by  deriving  values  of  k  from  the  spatially  atten¬ 
uated  spherical  waves  whose  propagation  is  described  by  the  Hufford 
equation.  Since  the  waves  from  which  we  construct  our  LORAN  pulse 
are  attenuated  (through  W) ,  they  have  a  spread  in  k  space.  That  spread 
introduces  some  uncertainty  into  the  coefficients  of  our  quadratic 
dispersion  relation.  The  uncertainty  in  turn  causes  a  spread  in  the 
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Table  1 

SPREAD  FACTORS  AND  PROPAGATION  VELOCITIES 
FOR  GROUND  CONDUCTIVITY  OF  0.01  mho/ra 
AT  SEVERAL  DISTANCES 


Distance 

Spreading 

Propagated 
rQ  (km) 

Factor 
(Ax/Ax0  -  1) 

V  -  V 
g  P 

(m/sec) 

76 

1.7 

X 

10-12 

4,779.3 

276 

0.9 

X 

10-10 

8,940.9 

476 

5 

X 

10-10 

11,656.6 

676 

1.5 

X 

io~9 

13,924.8 

876 

3 

X 

io-9 

16,356.6 

1,076 

6 

X 

io“9 

17,759.5 

2,076 

4 

X 

io’8 

25,155.6 

3,076 

1.4 

X 

io“7 

31,137.9 

NOTE:  The  group  velocity  is  Vg  and  the 
phase  velocity  is  Vp. 


Table  2 

SPREAD  FACTORS  AND  PROPAGATION  VELOCITIES 
FOR  GROUND  CONDUCTIVITY  OF  0.001  tnho/m 
AT  SEVERAL  DISTANCES 


Distance 
Propagated 
rQ  (km) 


Spreading 
Factor 
(Ax/Axq  - 


1) 


V  - 
g 


V 

P 


(m/sec) 


76 

2.5 

X 

10"9 

185,513.6 

276 

1 

X 

io’7 

-7 

315,607.6 

476 

3.4 

X 

10 

-7 

382,771.5 

676 

2.5 

X 

10 

389,607.4 

876 

1.5 

X 

io'12 

330,002.1 

1,076 

4.5 

X 

io-7 

-7 

249,938.0 

2,076 

2 

X 

10 

-8 

17,844.2 

3,076 

6 

X 

10 

2,897.4 

NOTE:  The  group  velocity  is  Vg  and  the 

phase  velocity  is  V  . 


computed  group  speed.  The  precision  to  which  we  are  attempting  to 
calculate  these  quantities  exceeds  the  precision  of  the  underlying 
concepts — phase  and  group  velocity. 

Unfortunately,  the  spread  in  group  velocity  for  lower  conductivity 
cases  can,  at  some  distances,  produce  uncertainties  of  the  order  of  the 
difference  in  the  phase  and  group  velocities.  Then  the  only  precise 
method  of  establishing  the  zero  crossing  times  is  to  transform  the 
LORAN  pulse  in  the  time  domain  to 'obtain  the  spectral  density,  multiply 
by  the  transfer  function  W,  and  finally  invert  the  transform  back  to  the 
time  domain.  The  last  step  is  known  to  be  too  complex  for  analytical 
solution.  However,  the  first  two  steps  of  the  procedure  produce  enough 
useful  information  to  be  worthwhile — e.g.,  how  the  frequency  band  in 
which  most  of  the  energy  in  the  pulse  is  carried  changes  with  distance 
traveled. 

Figures  4  through  6  illustrate  how  energy  partition  among  spectral 
components  changes  as  the  pulse  propagates.  Figure  4  gives  the  spectral 
distribution  at  the  transmitter  (derived  in  Appendix  B) .  Since  99  per¬ 
cent  of  the  energy  lies  between  90  and  110  kHz,  the  next  two  figures 
show,  for  ground  conductivities  of  0.01  and  0.001  mho/m,  how  the  energy 
balance  of  the  pulse  changes  between  76  and  3076  km,  respectively. 

The  attenuation  is  only  slightly  nonsymmetrical  about  100  kHz,  causing 
a  negligible  shift  toward  lower  frequency. 

CONCLUSIONS 

Based  on  numerical  computation  of  group  effects  on  a  LORAN  pulse, 
we  conclude 

•  A  low-conductivity  path  will  cause  a  significant  drift  of 
the  third  zero  crossing  within  the  pulse.  At  a  propagation 
distance  of  about  1300  to  2000  km,  cycle  ambiguity  can  enter. 

•  The  energy  remains  well  centered  at  100  kHz. 

Except  for  distant  receivers,  therefore,  the  time  of  arrival  of  the 
third  zero  crossing  is  given  accurately  by  the  prediction  of  the 


Relative  amplitude,  f(a>)  (dimensionless) 
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Hufford  equation  at  100  kHz.  "Distant11  depends  on  ground  properties, 
but  distances  up  to  about  1500  km  appear  to  present  little  difficulty, 
except  for  very  unusual  ground — e.g.,  Greenland  ice. 


Fig.  4--Spectral  distribution  of  transmitted  pulse 
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d  =  76  and  3076  km 
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Note:  This  figure  is  normalized  to  1 .0  at  100  kHz  for  both  distances.  Higher  frequencies 
are  attenuated  with  greater  distance.  The  gray  areas  show  the  relative  change  as  the 
pulse  propagates. 

Fig.  5--Amplitude  of  transfer  function  of  LORAN  pulse  at  76  and 
3076  km  for  ground  conductivity  0.01  mho/m 


Relative  amplitude,  f(a>)  (dimensionless) 
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d  =  76  and  3076  km 


Frequency  (kHz) 

Note:  This  figure  is  normalized  to  1.0  at  100  kHz  for  both  distances.  Higher  frequencies  are 
attenuated  with  greater  distance.  The  gray  areas  show  the  relative  change  as  the  pulse 
propagates.  This  figure  shows  greater  attenuation  than  Fig.  5. 


Fig.  6--Amplitude  of  transfer  function  of  LORAN  pulse  at  76  and 
3076  km  for  ground  conductivity  0.001  mho/m 
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IV.  NUMERICAL  STUDIES  OF  SURFACE  PHASE  PROPAGATION 


Here  we  discuss  the  results  of  the  PSR  computer  program  for  cal¬ 
culating  the  ground-wave  attenuation  function  W — which  describes  how 
the  ground  wave  signal  differs  from  free  space  propagation — according 
to 

ikr 

4>  =  2W  (const)  =  2W$0  ,  (11) 

where,  as  in  Sec.  Ill,  $  is  the  field  in  the  presence  of  the  ground 
and  is  the  free  space  solution. 

A  previous  PSR  report  [Field  and  Allen,  1978]  showed  that  VJ  for 
propagation  over  gentle  but  arbitrarily  variable  terrain  is  approx¬ 
imated  by  Huf ford's  Eq .  (11).  For  convenience,  we  repeat  that  equa¬ 
tion  here: 

mo  r°  r  -«l/2 

W(x0>  -  1  -  e-*’74  (i)  /  d*  [x-(- ~yj 

x  W(x) |  —  + 

\  n 
\  g 

where  =  the  distance  from  the  transmitter  to  the  receiver  pro¬ 
jected  on  a  flat  underlying  surface, 
k  =  the  free-space  wave  number  (jo/c, 
n  =  the  unit  vector  normal  to  the  ground, 

1/2 

n^  =  the  complex  index  of  refraction,  (tc  +  ia/wGQ)  , 
r^  =  the  straight-line  distance  between  the  transmitter  and 
the  receiver, 

=  the  straight-line  distance  between  the  transmitter  and 
the  integration  point  Q, 

r, ^  -  the  straight-line  distance  from  Q  to  the  receiver. 


3rA  ik(r1+r2-rQ) 


A  ik (r^4 

7  e 


Figure  7a  illustrates  the  geometry  of  r ,  r,,  and  r„. 
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a)  Hill  on  flat  earth 


z(x)  is  elevation 
measured  from  surface  S 


b)  Hill  on  smooth,  round  earth 


A  smooth,  round  earth  is  one 
Itrge  terrain  feature  dropping 
away  from  an  imaginary  flat 
surface  S  that  extends  from  the 
base  of  the  transmitter. 


Fig.  7--Geometrical  relationships  among  various  distances 
in  integral  equations  (12)  and  (13) 


Detailed  numerical  analysis  in  Field  and  Allen  demonstrated  that 
Eq.  (12)  gives  results  nearly  matching  the  well-known  residue  series 
for  a  smooth,  round  earth.  In  addition,  that  report  quantitatively 
examined  the  validity  criteria  for  Eq .  (12)  and  derived  a  more  exact 
integral  equation.  The  basis  of  the  PSR  ground-wave  computer  code, 
Eq.  (B-28)  of  Field  and  Allen  is  written 


W(xQ)  =  1  -  e 


x  W(x)  e 


1/2  r°  r  r2  11/2 

(*)  /  *  J 

0 
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The  improvements  over  the  Hufford  analysis  considered  here  are 
the  z"  impedance  correction  term  from  Sec.  II  and  the  following  dif¬ 
ferences  between  Eqs.  (12)  and  (13):  the  extra  ( i /kr?) (9r^ /8n)  term 

2  z  z 

in  Eq.  (13);  the  z*  (x)  term  in  Eq .  (13),  which  represents  an  arc- 

length  correction  over  undulating  ground,  zf  being  the  derivative  of 
the  ground  elevation  measured  with  respect  to  a  flat  surface;  and, 
finally  more  precise  use  of  the  distances  r^,  r^,  and  r^  instead  of 
their  projections  onto  a  flat  underlying  surface.  The  following  dis¬ 
cussion  presents  numerical  solutions  to  Eq .  (13)  for  generic  terrain 
characteristics  and  shows  how  inclusion  of  the  above  corrections 
changes  the  calculated  results. 

SMOOTH  EARTH 

The  simplest  generic  terrain  is  smooth — that  is,  without  local 
elevations  or  depressions  in  an  otherwise  flat  or  spherical  ground 
surface.  For  completeness  and  to  illustrate  basic  properties  of 
phase  delay,  we  first  examine  propagation  over  a  smooth  surface.  For 
a  smooth,  flat  surface,  z'  ,  and  z"  are  zero — as  are,  there¬ 

fore,  all  correction  terms.  To  show  how  the  phase  delay  of  100  kHz 
waves  depends  on  ground  conductivity  in  a  flat  earth.  Fig.  8,  calcu¬ 
lated  from  Eq.  (B-28)  of  Field  and  Allen,  displays  the  amplitude  and 
phase  of  W  for  conductivities  of  0.01  and  0.001  mho/m.  The  a  =  0.001 
case  gives  a  noticeably  stronger  effect. 

Figure  9,  also  calculated  from  Eq.  (B-28),  illustrates  the  same 
dependence  for  a  round  earth  with  no  correction  for  (i/kr^)  (S^/^n)  , 
zT,  or  zn.  These  three  functions  are  identically  zero  for  a  flat 
earth,  but  not  for  a  smooth,  round  earth — which  can  be  thought  of  as 
a  single  large  terrain  feature,  as  shown  in  Fig.  7b  above.  The  asso¬ 
ciated  corrections  are,  however,  too  small  to  show  significantly  on 
a  graph  for  a  smooth,  round  earth  with  a  =  0.01  mho/m. 

Figure  10  illustrates  the  difference  between  the  phase  of  W  for 
a  smooth,  flat  earth  and  that  for  a  smooth,  round  earth  with  all  of 
the  above  corrections.  Near  the  transmitter  the  difference  is  neg¬ 
ligible,  because  of  the  local  flatness  of  a  round  earth.  When  the 
earth’s  curvature  causes  the  ground  to  depart  by  an  amount  of  the 
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Fig.  8— Amplitude  and  phase  of  W  for  smooth,  flat  earth 
of  conductivities  0.01  and  0.001  mho/m 
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Fig.  9— Phase  of  W  for  smooth,  round  earth  of  conductivities  0.01 
and  0.001  mho/m  with  no  correction  for  (i/kr2)(3r2/3n) 
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Fig.  10— Phase  of  W  for  smooth,  round  earth  minus  phase 
of  W  for  smooth,  flat  earth 


order  of  a  wavelength  from  an  imaginary  flat  surface  extending  from 
the  base  of  the  transmitter,  we  expect  that  roundness  to  have  a  sig¬ 
nificant  effect.  For  a  100  kHz  signal  on  the  round  earth,  such  an 
effect  occurs  at  ~190  km.  The  ground  and  the  imaginary  surface  de¬ 
viate  by  a  few  percent  of  a  wavelength  at  ~35  km.  Thus  flat  and 
round  earth  results  should  vary  little  within  35  km  of  the  transmitter 
and  show  a  dramatic  difference  beyond  ~190  km. 

TERRAIN  ON  FLAT  EARTH 


Simple  Ramp  Terrain 


Terrain  influences  the  solution  to  the  integral  equation  (13) 
primarily  through  the  S^/Sn  term,  and  secondarily  through  the  cor¬ 
rection  terms  associated  with  Eq.  (13).  Therefore,  types  of  terrain 


possessing  simple  forms  of  9r2/9n  provide  useful  test  cases  from  which 
we  can  build  an  understanding  of  solutions  to  Eq .  (13)-  A  foil 
against  which  to  compare  our  numerical  results  for  simple  terrain  is 
what  we  hereafter  call  the  "f irst-order  solution" — which  is  valid  when 
terrain  is  sufficiently  gentle,  integration  paths  sufficiently  short, 
and  conductivity  sufficiently  high  that  W  may  be  considered  nearly 
constant  and  thus  can  be  factored  out  of  the  integral  on  the  right 
side  of  Eq.  (13).  The  integral  equation  can  then  be  solved  analytically. 

By  virtue  of  the  role  of  Sr^/^n  in  Eq .  (13),  the  effects  of  the 
terrain  features  are  most  readily  understood  in  terms  of  the  signs 
and  magnitude  of  the  first  and  second  derivatives  of  the  ground  ele¬ 
vation.  To  elucidate  the  important  effects  of  ground  shape  on  W,  we 
first  consider  terrain  features  that  are  a  sequence  of  flat  ramps, 
plotting  the  difference  between  the  phases  of  W  for  various  pairs  of 
surfaces — each  pair  comprising  a  particular  smooth  surface  and  the 
same  surface  with  a  terrain  feature.  We  hereafter  refer  to  that 
quantity  as  "phase  difference." 

Figure  11  plots  phase  difference  due  to  various  rises  of  constant 
slope  (ramps),  each  beginning  30  km  from  the  transmitter.  The  steeper 
slopes  exhibit  greater  phase  differences  at  the  foot  of  the  ramp  as 
well  as  greater  structure  in  certain  regions  along  the  ramp.  For 
moderately  steep  slopes,  that  property  may  be  explained  in  terms  of 
the  first-order  solution,  which  has  monotonic  behavior — slope  enter¬ 
ing  merely  as  a  multiplicative  factor.  For  much  steeper  slopes,  W 
changes  too  rapidly  for  the  first-order  solution  to  accurately  de¬ 
scribe  the  problem;  higher  order  terms  become  significant,  adding 
structure  to  the  ultimate  solution.  Near  the  foot  of  the  ramp — the 
point  at  which  the  second  derivative  of  the  ground  elevation  is  sig¬ 
nificant — the  phase  difference  is  strong. 

The  calculations  for  Fig.  11  ignore  the  (i/kr2) (Br^/Sn)  term. 

Its  importance  is  illustrated  for  a  ramp  of  slope  0.3  in  Fig.  12, 
which  plots  the  phase  difference  with  and  without  the  term.  For 
slopes  of  0.02  and  0.1,  its  effect  is  not  noticeable  on  a  plot  of 
this  scale.  Its  effect  is  strongest  near  slope  changes  and  weakest 
far  from  a  hill.  The  key  to  understanding  the  action  of  8r2/9n  is 
to  interpret  it  as  representing  the  angle  between  the  ground  at  an 
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Fig.  11— Phase  differences  for  several  constant  slopes  beginning 
at  50  km,  with  no  correction  for  (i/kr2) (3r2/3n) 


Fig.  12— Phase  difference  for  ramp  of  slope  0.3  beginning  at  50 
km,  with  and  without  correction  for  (1/kr2)(3r2/3n) 
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integration  point  and  a  straight  line  to  the  receiver.  When  the  inte¬ 
gration  point  and  receiver  are  on  the  same  flat  segment,  then 
is  zero.  On  flat  ground  beyond  a  region  of  changing  ground  slope, 
therefore,  S^/^n  is  zero  for  integration  points  everywhere  but  on  the 
rippled  region.  This  is  an  insignificant  part  of  the  path  except  for 
receivers  on  or  near  the  terrain  feature.  The  ( i /k^) (Br^/Bn)  term 
introduces  small-scale  structure  near  regions  of  strong  slope  change. 
For  example,  numerical  calculation  shows  that  at  the  foot  of  a  0.1 
upward  slope,  a  phase  difference  equivalent  to  85  m  appears  in  tin 
first  0.1  km,  which  would  be  missed  either  by  ignoring  this  term  or 
by  calculating  intervals  too  wide  to  focus  on  the  region  of  large 
second  derivative.  The  ( i/kr^) (Br^/Sn)  term  is  included  in  all  sub¬ 
sequent  figures  unless  otherwise  specified. 

The  stronger  effect  on  phase  of  a  downward-slop i ng  ramp  relative 
to  that  of  a  ramp  of  equal  and  opposite  slope,  shown  in  Figs.  13  and 
14,  cannot  be  explained  by  the  first-order  solution  since,  at  that 
order,  the  phase  difference  is  proportional  to  the  slope  of  the  ramp. 
The  stronger  effect  is  due  to  the  magnitude  of  W  working  together  with 
in  the  integral  to  produce  the  actual  solution.  Figure  15 
displays  the  amplitude  of  W  associated  with  plus  and  minus  0.3  slope 
ramps  that  begin  50  km  from  the  transmitter.  It  shows  the  source  of 
the  phase  difference  structure  exhibited  in  Fig.  11  for  0.3  slope. 

The  plot  of  the  amplitude  of  W  for  the  downward-sloping  ramp  also  ex¬ 
hibits  the  monotonic  phase  difference  structure  evident  in  Fig.  13 
for  -0.3  slope. 

The  arc-length  correction  for  a  0.1  slope  ramp  is  about  1  m  in 

apparent  propagation  distance,  and  is  thus  too  small  to  show  clearly 
* 

on  the  graphs.  On  more  complex  terrain,  the  arc  length  may  have 
greater  effect.  Its  influence  on  phase  is  less  than  that  of  stretch¬ 
ing  the  propagation  path  out  flat  and  including  the  extra  distance 
traveled,  since  the  phase  of  W  is  a  nonlinear  function  of  distance 
obtained  by  taking  the  arc-tangent  of  the  ratio  of  Im(W)  and  Re(W) . 

*The  plots  of  phase  difference  scales  of  0.2tt  can  resolve  phase 
differences  of  about  0.002  rad.  Since  477  times  the  phase  difference 
in  radians  gives  the  distance  difference  in  meters,  one  meter  would 
appear  to  merely  thicken  linos  on  the  plot. 
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Fig.  13--Phase  difference  for  ramps  of  slopes  -0.1  and  -0.3 
both  beginning  at  50  km 


Fig.  14--Phase  difference  for  ramps  of  slopes  -0.1  and  0.1 
both  beginning  at  50  km 
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Fig,  15--Amplitude  of  W  for  ramps  of  slopes  0.3  and  -0.3 


We  next  consider  two  examples  showing  what  happens  when  the  slope 
of  the  ground  changes  abruptly.  Our  calculations  hereafter  include  the 
arc-length  correction.  Figure  16  plots  the  phase  difference  asso¬ 
ciated  with  a  ramp  of  slope  0.1  beginning  at  50  km  and  changing  to 
a  0.2  slope  at  100  km.  The  slope  of  the  phase  difference  does  not 
change  sign  at  the  junction  of  the  two  slopes.  Between  50  and  100  km, 
the  difference  matches  the  behavior  of  the  0.1  ramp  found  previously 
in  Fig.  11.  After  100  km,  the  phase  difference  resembles  that  found 
for  the  0.2  slope  in  Fig.  11,  but  its  structure  is  somewhat  more 
pronounced . 

The  behavior  of  the  second  case  is  strikingly  different.  Fig¬ 
ure  17  plots  the  phase  difference  associated  with,  again,  a  0.1  ramp 
beginning  at  50  km  but  changing  to  a  0.05  slope  at  100  km.  The 
second  slope  is  positive  and  still  represents  an  upward  ramp,  but 


difference  i  Phase  difference 


j.  16--Phase  difference  for  ramp  beginning  at  50  km  that 
changes  slope  from  0.1  to  0.2  at  100  km 


Fig.  17--Phase  difference  for  ramp  beginning  at  50  km  that 
changes  slope  from  0.1  to  0.05  at  100  km 
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not  so  steep  as  the  immediately  preceding  0.1  slope.  Indeed,  the 
derivative  of  the  phase  difference  plot  changes  sign  at  the  junction 
of  the  two  slopes.  Whereas  the  phase  difference  plot  behaves  in  the 
established  pattern  up  to  100  km,  the  difference  associated  with  the 
0.05  slope  looks  not  like  that  of  an  upward  ramp,  but  like  that  of  a 
downward  ramp — which  it  is  relative  to  the  preceding  slope .  The 
integral  equation  thus  exhibits  a  "memory"  of  the  ground  terrain,  so 
that  the  change  in  slope  is  a  governing  feature  of  how  phase  differ¬ 
ences  arise. 

1/2 

Referring  back  to  Eq.  (13),  we  see  that  the  presence  of  (1/^) 
in  the  integrand  tends  to  make  the  slope  immediately  preceding  the 
receiver  most  important,  whereas  slopes  farther  from  the  receiver  con¬ 
tribute  weakly.  That  sensitivity  to  slope  change  can  be  regarded  then 
as  sensitivity  to  the  sign  of  the  second  derivative  of  the  ground  ele¬ 
vation.  Thus  we  see  a  consistency  in  the  extra  structure  of  the  phase 
difference  for  the  0.2  slope  part  of  Fig.  16  and  the  inverted  behavior 
of  the  0.05  slope  part  of  Fig.  17. 

Finally,  considering  a  one-dimensional  pyramid  having  equal  ascend 
ing  and  descending  slopes,  Fig.  18  plots  the  characteristic  phase  dif¬ 
ference  pattern  for  a  symmetrical  hill — i.e.,  a  phase  decrease  (over 
the  rising  side)  followed  by  a  sharp  increase  (near  the  peak  on  the 
down  side)  followed  by  a  long  tail.  Similar  patterns  have  been  found 
by  others  (e.g.,  Johler  and  Berry  [1967])  for  different  symmetrical 
one-dimensional  hills.  The  existence  of  the  long  tail  is  significant 
for  navigation,  since  it  indicates  that  a  hill  not  only  causes  a  local 
effect,  but  that  large  phase  differences  (perhaps  many  tens  of  meters, 
depending  on  the  size  of  the  hill)  persist  many  hill  lengths  beyond 
the  hill  itself.  The  tail  occurs  because  the  second-derivative  effect 
just  discussed  prevents  the  effects  of  the  ascending  and  descending 
sides  of  the  hill  from  cancelling. 
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Fig.  18--Phase  difference  for  pyramid  beginning  at  50  km  with 
faces  of  slopes  0.1  and  -0.1  (characteristic 
signature  of  symmetrical  hill) 


slope.  We  use  a  1  *  10  sin  (x)  elevation  function  that  joins  smoothly 

k 

with  the  flat  ground.  The  second  derivative  changes  sign  on  each 
side  of  a  given  hill.  Figure  19  shows  the  results  of  including  the 

2 

( i  /k^)  (Br^/Sn)  ,  arc-length,  and  z"  impedance  corrections  for  a  sin 
hill  of  roughly  the  same  proportions  as  the  pyramid  in  Fig,  18,  Apart 
from  somewhat  less  smooth  transitions  at  certain  places,  the  pyramid 
behaves  much  the  same  as  the  more  complicated  sin^  hill. 

Figure  20  directly  contrasts  the  effect  of  the  same  two  hills  on 
phase.  It  is  notable  that  merely  matching  height  and  width  of  the 
two  hills  yields  a  close  match  in  the  traces  of  the  phase  differences. 


k 

We  describe  hills  by  giving  the  height  times  the  half-width  in 
kilometers,  unless  otherwise  specified.  We  abbreviate  the  height 
function  of  a  hill  7.(x)  =  il  sin2(x  -  x0/L)  as  II  *  L  sin2.  Thus  a 
,fl  *  10  sin  "  lull  has  n  peak  1  km  high,  a  half-width  of  10  km,  and 
a  sin^  profile. 


Phase  difference 
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Fig.  19--Phase  difference  for  1  x  10  sin  hill  beginning  at 
50  km,  with  and  without  correction  for  (i/kr2) (3r2/3n) 


Fig.  20--Phase  difference  plot  comparing  1  x  10  pyramid 
with  1  x  10  sin^  hill,  both  beginning  at  50  km 


Figure  21  illustrates  how  the  phase  difference  changes  over  the  hill 
when  the  conductivity  diminishes  by  a  factor  of  10  to  0.001  mho/m. 
Although  1/n  is  much  larger  than  before,  Sr^/Sn  is  still  the  dominant 
term  near  the  hill,  so  the  only  observable  difference  is  in  the  tail 
and  at  the  negative  phase  peak. 

TERRAIN  ON  ROUND  EARTH 

Hills 

Effect  of  Shape  of  Hill.  Having  explored  the  properties  of  hills 

alone,  we  couple  the  hills  with  round-earth  propagation.  Figure  22 

2 

shows  three  different  sin  hills  50  km  from  the  transmitter  on  a  round 

2 

earth.  The  effect  of  the  arc-length  correction  on  the  1  x  10  sin  hill 
is  about  2  m  in  apparent  propagation  distance,  and  thus  too  small  to 
show  clearly  on  the  graphs. 

2 

The  z"  impedance  correction  is  negligible  even  for  a  1  *  3  sin 
hill.  We  may  safely  ignore  this  factor  because  the  impedance  correc¬ 
tion  will  be  of  consequence  only  when  the  radius  of  a  hill  becomes 
comparable  to  the  skin  depth  (tens  of  meters);  moreover,  no  reliable 
topographical  data  of  such  precision  are  available. 

Effect  of  Distance  to  Hill.  Now  we  explore  the  effect  of  hill 
distance  from  the  transmitter,  considering,  first,  hills  very  close 
to  the  transmitter,  within  the  distance  where  earth  curvature  effects 
should  be  negligible;  second,  those  at  the  horizon;  and  third,  those 
clearly  over  the  horizon.  A  hill  is  considered  on  the  horizon  when 
the  line  of  sight  from  the  transmitter  tangent  to  the  earth  touches 
the  top  of  the  hill.  At  110  km,  a  hill  1  km  high  sees  the  transmitter 
as  being  just  on  the  horizon. 

Figure  23  illustrates  the  subtle  evolution  of  the  signature  of  a 
symmetrical  hill  with  increasing  distance  from  a  transmitter  on  a 
round  earth.  Most  noticeably,  the  tail  is  larger  when  the  distance 
to  the  transmitter  is  smaller.  The  same  effect  would  appear  on  flat 
ground  to  a  similar  degree.  For  closer  hills,  the  initial  phase  de¬ 
crease  is  shallower  and  the  peak  leading  into  1 1  le  tail  is  higher. 

1/2 

This  behavior  is  due  primarily  to  the  (l/r^)  term  in  the  inte¬ 
grand  of  Eq .  (13),  which  acts  as  a  distance  weighting  factor  and 
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Fig.  21--Phase  difference  for  1  x  lo  sin  hill  beginning 
at  50  km,  with  conductivities  0.01  and  0.001  mho/m 


Fig.  22 — Phase  differences  for  three  different  sin^  hills 
beginning  at  50  km 
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Fig.  23— Phase  differences  for  1  x  10  sin2  hills  beginning 
at  5,  110,  and  400  km 


causes  hills  near  the  transmitter  to  have  a  stronger  tail  effect  than 
those  farther  away. 

We  do  not  show  the  analogous  plot  for  the  same  hills  on  a  flat 
earth  because  it  Is  nearly  identical  to  Fig.  23.  There  is  a  small 
difference  in  the  behavior  of  hills  on  round  and  flat  earth  because, 
on  round  earth,  the  slope  of  a  face  of  a  hill,  relative  to  horizontal 
at  the  transmitter,  depends  on  the  slope's  distance  from  the  trans- 
mitter,  as  is  illustrated  in  Fig.  24. 

Proximity  to  the  horizon  seems  to  have  no  particular  effect, 
because  the  receiver  is  on  the  ground.  That  lack  of  effect  may  be 
readily  understood  in  terms  of  a  division  of  the  wave  field  into  space 
waves  and  surface  waves.  Norton  [1937]  divided  the  field  into  a  part 
that  drives  the  ground  currents  and  one  that  does  not,  and  found  that 
the  first  part  drops  off  with  height  above  the  surface.  That  part  is 
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Fig.  24--Geometry  of  hills  near  and  far  from  transmitter 


the  surface  wave.  The  other  part  decreases  to  zero  at  the  surface 
and  increases  with  height.  The  horizon  has  no  effect  because  the 
surface  wave  component  so  strongly  dominates  the  ground  wave  at  small 
heights  that  loss  of  the  direct  wave  component  upon  passing  the  horizon 
is  unimportant. 

Valleys 

We  next  consider  the  effect  on  phase  of  a  valley  that  is  on  the 

2 

propagation  path.  Figure  25  compares  the  effect  of  a  1  *  10  sin  hill 

2 

with  that  of  its  -sin  complement — i.e.,  valley.  Although  the  general 
shapes  of  the  phase  difference  plots  are  reflections  about  zero,  both 
hill  and  valley  have  a  positive  tail.  For  a  moderately  sloped  valley, 
such  as  that  in  Fig.  25,  the  incomplete  mirroring  of  the  phase  dif¬ 
ference  of  a  hill  by  that  of  Us  complementary  valley  can  be  explained 
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Fig.  25--Phase  difference  plot  comparing  1  *  10  sin  hill 
and  valley,  both  beginning  at  50  km 


1 

\ 

i 

in  terms  of  our  qualitative  slope  and  second-derivative  arguments.  • 

2  ! 
The  ground  elevation  function  of  a  sin  hill  has  two  regions  of  posi¬ 
tive  second  derivative  and  only  one  of  negative  second  derivative; 
the  valley,  however,  has  two  regions  of  negative  second  derivative 
and  only  one  of  positive  second  derivative.  As  discussed  earlier, 

moderate  negative  second  derivatives  more  strongly  influence  the  phase  j 

difference  than  equal  positive  second  derivatives.  Such  lack  of 
symmetry  in  the  effects  of  a  moderate  valley  and  the  corresponding 

moderate  hill  accounts  for  the  qualitative  discrepancy  in  phase  dif-  , 

ference. 

Slightly  and  steeply  sloped  valleys  have  different  effects.  For 

2 

example,  a  weak  1  x  50  sin  valley  produces  a  phase  difference  that 
corresponds  qualitatively  to  the  predictions  of  a  first-order 


i 

j 

! 


calculation  and  has  a  negative  tail,  as  shown  in  Fig.  26.  At  the 

opposite  extreme,  we  consider  a  hill  and  valley  given  by  a  7.5  x  25 
2 

sin  ground  elevation  function,  shown  in  Fig.  27.  These  features 
correspond  to  Mt .  Everest  and  a  valley  of  the  same  size.  The  steep 
slope,  long  faces,  and  attendant  large  second  derivatives  combine  to 
give  phase  differences  somewhat  distorted  from  the  pattern  we  have 
found  thus  far.  This  last  case  also  stretches  the  limit  of  applic¬ 
ability  of  the  integral  equation  because  the  approximation  of  a  Hertz 
potential  with  a  single  component  normal  to  the  earth  is  not  very 
accurate  over  large  portions  of  the  integration  path  [Field  and  Allen, 
1978]. 

Hill  and  Valley  Combinations 

Next  we  consider  combinations  of  a  hill  and  a  valley.  Figures  28 
through  31  give  the  characteristic  signatures  of  a  hill  followed  immed¬ 
iately  by  a  valley,  a  valley  followed  by  a  hill,  and  the  same  combina¬ 
tions  with  a  smooth  region  the  size  of  the  hill  separating  the  two. 
Rather  than  any  tendency  for  each  to  cancel  the  other's  effects,  the 
plots  indicate  that  the  valley  and  hill  keep  the  integrity  of  their 
signature  when  isolated,  merely  beginning  on  the  tail  of  the  preceding 
terrain  feature.  Notably,  the  major  effect  of  a  hill/valley  combina¬ 
tion  is  not  merely  the  extra  distance  traveled  by  the  wave.  This  is 
shown  by  computing  the  phase  delay  with  and  without  the  arc-length 
correction.  The  difference  accounts  for  approximately  2  m  in  the 
distance  determined  from  phase  difference. 


Minimum  Significant  Terrain  Feature 

Next  we  look  for  the  minimum  terrain  feature  that  has  a  signifi¬ 
cant  effect,  in  order  to  determine  what  data  filtering  is  acceptable 
before  applying  a  one-dimensional  integration  technique.  A  small 
individual  terrain  feature  has  a  small  effect.  Because  of  the  tendency 
of  features  to  build  onto  the  tail  of  the  previous  feature,  however, 
long  chains  of  small  hills  can  result  in  large  accumulated  phase  devia¬ 
tions,  even  if  the  hills'  height  and  width  are  in  the  same  1:10  ratio 
we  have  been  studying.  Moreover,  even  for  hills  only  100  m  high  and 
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Fig.  28— Phase  difference  for  lx  10  sin2  hill  beginning  at 
50  km,  followed  by  1  x  10  sin2  valley  beginning  at  70  km 


Fig.  29— Phase  difference  for  1  x  lo  sin2  valley  beginning  at 
50  km,  followed  by  1  x  10  sin2  hill  beginning  at  70  km 
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Fig.  30--Phase  difference  for  1  x  10  sin  hill  beginning  at 
50  km,  followed  by  1  x  10  sin2  valley  beginning  at  90  km 


Fig-  3,1-p^a?®  difference  for  1  x  10  sin2  valley  beginning  at 
50  km,  followed  by  1  x  10  sin^  hill  beginning  at  90  km 
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1000  m  long,  a  sequence  of  10  hills  can  produce  a  maximum  phase  dif¬ 
ference  equivalent  to  approximately  40  m.  Longer  sequences  produce 
larger  accumulated  phase  differences.  Figures  32  through  34  illustrate 
such  behavior. 

The  present  analysis  indicates  that  a  sufficiently  long  chain  of 
hills — however  small  the  individual  hills — will  produce  a  noticeable 
phase  shift.  This  result  is  not  to  be  expected  of  real  terrain,  which 
generally  has  transverse  and  longitudinal  extensions  of  similar  di¬ 
mensions.  Isolated  ridges  are  common,  but  not  extended  series  of 
ridges.  The  small  features  thus  are  not  realistically  represented  by 
the  present  one-dimensional  treatment.  Field  and  Joiner  [1979]  have 
shown  that  the  effect  of  an  obstacle  occluding  a  small  portion  of  a 
Fresnel  zone  is  reduced  pro  rata  from  that  of  an  obstacle  filling  the 
Fresnel  zone.  Extension  to  infinity  in  the  transverse  direction  has 
little  more  effect  than  just  filling  the  Fresnel  zone. 

More  explicitly,  the  one-dimensional  approximation  (13)  to  the 
full-fledged  two-dimensional  equation  (C.l)  assumes  no  ground  structure 
in  the  transverse  direction.  The  phase  difference  for  hill  chains 
just  calculated  is  analogous  to  propagation  across  the  ridges  of 
washboard- 1 ike  terrain  wide  enough  to  fill  a  Fresnel  zone  of  width 

/(wave  length) x  (receiver  distance  from  transmitter)  . 

The  width  of  a  Fresnel  zone  at  200  km,  for  example,  is  v3  km  y  200  km 
25  km.  A  chain  of  hills  1  km  long  and  25  km  wide  is  unlikely  to 
extend  many  kilometers.  One  would  be  more  likely  to  find  terrain 
with  many  scattered  hills  about  1  km  long  and  about  1  km  wide.  A 
single  such  hill  would  have  roughly  1/25  the  effect  of  a  hill  1  km 
long  and  25  km  wide.  Of  course,  the  hill  25  km  wide  would  fill  the 
Fresnel  zone  and  have  essentially  the  same  effect  as  an  infinitely 
wide  hill.  More  precise  evaluation  of  the  effect  of  obstacles  small 
relative  to  the  width  of  a  Fresnel  zone  requires  recourse  to  the  two- 
dimensional  integral  equation  (C.l). 

CONCLUSIONS 

We  have  presented  numerical  results  for  a  variety  of  possible 


A 
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Fig.  32--Phase  difference  for  sequence  of  two 
0.5  x  5  sin?  hills 


Fig.  33--Phase  difference  for  sequence  of  four 
0.25  x  2.5  sin2  hills 


-43- 


Fig.  34--Phase  difference  for  sequence  of  ten 
0.1  x  1.0  sin2  hills 

terrain  configurations  with  an  eye  toward  determining  which  factors 
are  significant  in  modifying  the  calculated  time  of  arrival  of  phase 
for  a  100  kHz  signal.  We  conclude  the  following: 

•  Terrain  features  produce  local  and  distant  perturbations 

in  propagation  of  phase  relative  to  a  smooth  surface.  Such 
effects  are  characteristic  of  the  type  of  terrain  feature 
encountered  (e.g.,  hill,  valley)  and  can  be  equivalent  to 
many  tens  of  meters  of  recei ver-to- transmi t ter  distance  on 
smooth  ground  even  for  gentle  terrain  over  distances  of  a 
few  hundred  kilometers. 

•  The  significance  of  arc  length  and  curvature-modified  im¬ 
pedance  corrections  is  small,  amounting  to  a  few  meters  under 
typical  conditions.  The  inclusion  of  the  ( i/k^)  (S^/Sn) 
term  is  important,  especially  in  regions  of  large  slope 
change  such  as  at  the  base  of  a  hill. 
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•  The  magnitude  of  the  effect  of  a  particular  terrain  feature 
is  insignificantly  affected  by  whether  the  underlying  smooth 
surface  is  flat  or  round. 

2 

•  Comparing  Fig.  19,  which  illustrates  a  sin  hill  and  a 
pyramid  of  the  same  length  and  height,  with  Fig.  9,  which 
shows  the  effect  of  different  ground  conductivities  on  phase, 
we  find  that  knowledge  of  ground  conductivity  is  more  critical 
than  knowledge  of  the  exact  shape  of  a  terrain  feature. 

•  Sequences  of  terrain  features  yield  phase  differences  roughly 
equivalent  to  those  resulting  from  starting  the  phase  dif¬ 
ference  patterns  of  subsequent  bumps  on  the  tail  of  preceding 
ground  features.  Since  such  an  effect  is  cumulative,  many 
small  bumps  will  significantly  alter  the  effective  distance 
as  determined  by  phase  arrival.  The  one-dimensional  integral 
equation  (13)  will  overestimate  the  effect  of  such  terrain, 
which  will  rarely  fill  the  entire  width  of  a  Fresnel  zone. 
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Appendix  A 


DERIVATION  OF  CURVATURE  CORRECTION  TO  SURFACE  IMPEDANCE 


Performing  a  detailed  derivation  of  Eq.  (8)  (Sec.  II,  p.  5),  we 
calculate  here  the  impedance  e/H  for  the  case  of  a  wave  propagating 
at  the  surface  of  a  curved  conductor,  with  a  correction  for  curvature. 
Figure  1  (p.  5)  illustrates  the  geometry  of  the  problem. 

We  proceed  from  the  following  two  Maxwell  equations: 

-V 

V  X  E  =  -v  |~  ,  (A. 1) 


and 


V  x  H  =  at 


(A. 2) 


We  assume  the  fields  to  be  of  the  following  forms,  respectively: 


E  =  ea)y  exp 


-  <E0+  SEj  +  «2E2  ♦  ,  (A.3) 


and 


H  =  h  exp  <  i 


2)-  -  u)t]| 


=  (Hq  +  6H1  +  62H2  + 


(A. A) 


The  skin  depth  6  is  given  by 


(A. 5) 


PRECEDING  PAGE  BLANK-MOT  FILMED 


Equation  (A.l)  thus  becomes 


which  may  be  written  as 


V  x  e  - 
x 


V  x  e  - 

y 


V  X  e  - 
z 


Similarly  rewriting  Eq . 


V  x  h  = 
x 


V  x  h  = 

y 


V  x  h  = 
z 


V  x  E  =  iojpli  , 

independent  component  equations: 
I  <ey***'  '  e^’y,)  =  ihx  ’ 

i  (<V*S'  ~  =  ihy  * 

I  (ex*V  "  =  lhz  • 

(A. 2)  yields 

I  <V’z’  "  h**’y')  =  Wyex  ’ 

i  (V’x'  '  hx<V>  =  Wayey  ’ 

t  <hX^’y-  -  V’x,:>  =  • 


(A. 6) 


(A.  7) 


A  prime  on  a  differentiating  subscript  denotes  that  the  coordinate 

is  in  units  of  skin  depth,  e.g.,  y*  =  y/6. 

We  now  describe  the  surface  locally  by  two  unit  tangent  vectors, 
-y  ~y  -> 

A  and  B,  and  a  normal  n  pointing  into  the  conductor  with  components 

(n  ,  n  ,  n  ).  The  condition  of  normal  propagation  is  that  ^  be 
x  y  z  _  ^ 

along  n  so  that  ^  x  n  =  0.  Therefore,  we  now  write 


n 

ijj,  = 

y  7.  n 


For  simplicity  of  notation,  hereafter  we  use  a  comma  followed 
by  a  subscript  to  denote  partial  differentiation,  e.g.,  =  9^/3z. 
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ip,  =  ip,  ~ 

x  z  n 


From  Eq.  (A. 5)  it  follows  that  o  *  2/(iop6  ),  so  we  may  rewrite 
with  an  explicit  ordering  of  6  the  equation  sets  (A. 6)  and  (A. 7): 

(n  . 

e  i|»,  t  -  e, ip,  ,  -Z )  *  i6h  ,  (A. 8a) 

yr  *  z  zY  z  n  J  x 


6(V  x  e) 

y 


-  1  (ez^z-  f  '  eJ’z')  =  i<5hy  ’ 


6(Vz  x 


-  1  (ex*»*'  5*"  ey^z-  ^)=  *«*,  . 
z  z  ' 


«2«, *  s>  - iS  ( %♦  V  -  /)  ■  2%  • 

‘!'V*)  - 16  • 

/  n  n  . 

62 (V  x  h)  -  ifi  h  4),  ,  -  h  \p,  ,  —  1  =  2e 

x  z  \  xY  z  n  y  z  n  / 


6  (V  x  h) 

y 


The  6  terms  yield 


■*  (  VKz'  ■  ‘d,*'.'  /)  "  0  • 

-*  ( V-z'  r  -  Vs-)  -°  • 

'  z 

l  n  n  \ 

-i  (  E.  tf>,  ,  -  E.  {(1,  ,  —  I  =  0  . 

\  0xy  z  n  Oy  z  n  ) 


(A. 8b) 


(A. 8c) 


(A. 9a) 


(A. 9b) 


(A. 9c) 
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For  a  hill,  as  shown  in  Fig.  1  (p.  5) ,  n  =  n  =  n  is  excluded. 

x  y  z 

The  last  set  of  equations  then  implies 


E0x  E0y 


E0z  =  0 


(A. 10) 


Proceeding  to  the  terms  of  order  6^  that  appear  in  Eqs.  (A.8a-c) 

■+  —v 

A.9a-c),  it  is  found  that  n  • 
terms  may  also  be  combined  to  obtain 


and  (A.9a-c),  it  is  found  that  n  •  ==  0  and  n  #  E^  =  0.  The  same 


n  «  n  n  n 

H0z  r  *’  +  h 


n  n 

2  r  “Oy 


n 


Ox 


r2  [  i  +  i  -  2i 


Oz 


f  n 

_  JL 

n 


,,2' 


+  H 


0y 


-21  +  ip 


,2 


1  + 


x 
2 


Ox 


n  n 

2 

n 


2 


(A. 11) 


where  If  we  require  these  equations  to  bo  consistent  for 

all  values  of  and  H„  ,  then  the  respective  coefficients  of 

(Jy  Oz  0v 

and  H„  must  be  equal.  That  equality  leads  to  the  condition 
Oz 


n 

7. 


(A. 12) 
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But 


^  •  n  =  *,x,nx  +  *,y,ny  +  *z 


n2  n2  n2 

=  f  +  /  +  n* 

z  z  z 


.  41 

n 

2 

The  derivative  of  the  phase  \Jj  in  the  normal  direction  is  thus  ±/2i . 
Considering  the  refractive  material  to  be  a  hill  with  free  space 
above,  we  want  the  field  to  attenuate  with  decreasing  z.  We  thus 
choose  the  4*  sign  in  Eq .  (A.  12). 

Because  of  Eq.  (A. 12),  the  field  resulting  from  terms  of  order 
6  obeys 


J1A 


/2i 


H0B  ~  Hqb 


(A. 13a) 


and 


E1B  ■  V  '  “OA  •  <A-nb) 

These  two  equations  are  obtained  simply  by  projecting  the  6^  result 
- ►  -> 

onto  the  A  and  B  basis  vectors. 

->  ->  2 

Then,  taking  projections  onto  the  A,  B  basis  of  the  0  terms  and 
using  previous  results  to  simplify,  after  lengthy  calculation,  we 
obtain  the  following  set  of  equations: 
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(V  x 

(V  x 

(V  x 

(V  x 

(V  x 

(V  x 


h\ 


V  B 


Vn 


Va 


Vb 


H„ ) 
On 


. 

-  1  r  E2B  -  1H1A  ’ 
2 


.  V 

+  1  n  -  E2A  =  1H1B  * 

Z 


"  iHln  * 


.  V 

-  1  n  "iB  ■  2E2A  ■ 

z 


+  i  ^  H  =  2E  , 
n  1A  2B 
z 


=  2E 


2n  • 


(A. 14a) 

(A. 14b) 

(A. 14c) 

(A. 14d) 

(A. 14e) 

(A. 14f ) 


If  Eq .  (A. 14)  is  to  have  a  self-consistent  solution,  then  (A. 14a)  and 

(A.14e)  must  yield  (V  x  E.)  =  i/i/2  (V  x  H  )  ,  (A. 14b)  and  (A.14d) 

1  A  U  d 

must  yield 


(V  X  E1)B  =  v'iTI  (V  X  H0) 


A  ’ 


(A. 15) 


and  (A.14d)  must  yield  E0  =  -i/i/2  H,_  +  1/2(V  x  H  )  . 

^A  Id  U  A 

Recall  that  the  electric  field  of  our  wave  is 


EA  ■  SE1A  +  J  E2A  +  - 


Substituting  for  and  E^^  from  Eqs.  (A. 13a)  and  (A. 13), 


Ea  -  6(-i/i?2  H0fi)  +  62[-i/i72  Hln  +  1/2 (V  x  ha)a  +  •••]  .  (A. 16) 


IB 


0  A 


The  H^b  term  is  associated  with  the  curvature  of  the  incident  wave 

front  and,  for  this  problem,  we  let  H,_  =  0.  We  evaluate  1/2(V  *  iL). , 

Id  U  A 

assuming 
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Appendix  B 

CALCULATION  OF  LORAN  SPECTRAL  INTENSITY  FUNCTION 


To  analyze  group  effects,  we  must  know  the  intensity  at  all  fre¬ 
quencies  of  a  LORAN  pulse.  Here  we  derive  the  functional  form  of  the 
LORAN  pulse  in  the  frequency  domain  plotted  by  Fig.  4  (Sec.  Ill,  p. 

15) — which,  when  multiplied  by  W,  is  the  basis  for  Figs.  5  and  6. 

A  single  LORAN  pulse  has  the  form 

/  ,  \2  2 ( 1—  t / 1  ) 

F ( c )  =  I  — )  e  p  sin  (a^t  +  0)  ,  x  =  0,  t  a  0  .  (B.l) 

V  p/ 

We  assume  that  each  pulse  effectively  decays  to  zero  before  the  emis¬ 
sion  of  the  following  pulse.  We  further  assume  that  propagation 
effects  are  not  so  severe  as  to  ever  cause  pulses  to  overlap.  Thus 
we  may  analyze  only  a  single  pulse.  The  spectral  intensity  is 


f  (u),  x  -  0)  = 


CO 

(27 r)1/2  / 


.  -iu)t  , 
F(t)e  dt 


(B .  2) 


If  we  define  step(t)  =  1  for  t  >  0,  0  otherwise,  then 


f(w.  0)  -  - 


1  f  T  t  (1"t/tn)l" 

/ step<t)  [-, 0  J 


sin  (u^t  +  0)e  dt 


.  2  f  7  -2t/t 

- ~~k  I  step ( t )  t  e  p  sin  (u)nt  +  0)e  1U)t  dt 

(2tt)  l' z  t  J  u 

v  p 


c* 

_i_- .  £-2  f 

(277)  1/2  t?  / 


—  (  iaH-2  / 1  )t 


t  sin  (to  t  +  0)e 


P  0 


dt  . 


(B.3) 


Then 


1  ^  n 

f(0),  0)  =  - ~nr^£[t  sin  (ft)  t  +  0)) 

(2  ft))1'  t 

P 


=  - ^  ^  ^-Xtsin  (0)  t  +  0)]  .  (B.6) 

(2ir)x/^  t  dp 
P 


Using 


sin  (od^t  +  0)  =  sin  co^t  cos  0  +  cos  u^t  sin  9  ,  (B.7) 


we  obtain 
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Now 


(B.13) 


1 
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Let  us  take  0=0.  Then 


f(w.  0)  = 


_  2p  _  1  e2  /  2  2  \ 

(>2  +  Sy& 7/1?r 


Using 


P  =  iw  +  —  > 
P 


2  2  i4co  4 

p  .  -a.  +  —  +  -2  • 

P  t 

P 


we  obtain 


Ref  (co 


if),  .(is)*] 


and 


Imf  (u>,  0)  -  n  j^~  ro  +  ^-3w2  -  a>2  +  A| 
1  ‘  1P 


)*]• 


where 


2  2 
2e  / t 

P 


■t[0 

p  p '  '  '  J 


and 


(B.14) 


(B.15) 


(B. 16) 


i 
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Appendix  C 

CONSISTENCY  OF  HEIGHT  GAIN  FUNCTION  FORMULATIONS 
AT  FLAT  PLANE  SURFACE 


The  attenuation  function  for  waves  propagating  on  a  rippled  sur¬ 
face  is  given  by  Hufford  [ 19  52  j  to  satisfy 


W(P) 


+  ik- 

2ir 


f* 


_o 

'lr2 


W(Q) 


ik(r1+r2~r0) 


(C.l) 


When  the  receiver  is  above  the  surface,  however,  the  Hufford  treat¬ 
ment  gives 


2W(P) 


rir2 


W(Q) 


ik(r1+r2-r0) 


\  9l2l 

J 


(C.2) 


Equations  (C.l)  and  (C.2)  appear  to  be  identical,  except  for  a 
factor  of  2  in  (C.2).  Therefore,  one  would  expect  that  at  t h*:  sur¬ 
face  of  the  earth — i.e.,  when  the  elevation  above  the  earth  is  reduced 
to  zero — the  solution  of  W  would  depend  on  which  equation  one  chose  to 
use.  Because  physical  phenomena  tend  to  behave  smoothly,  however,  a 
1  m  elevation  should  not  significantly  change  a  system  with  a  3000  m 
wavelength.  Therefore,  we  look  more  closely  at  Eqs.  (C.l)  and  (C.2). 

This  appendix  demonstrates  mathematically  that  the  two  equations 
have  additional  differences,  so  they  actually  give  a  continuous  answer 
at  the  surface  of  the  earth.  Two  of  the  differences  are  important. 
First,  Eq .  (C.l)  is  an  integral  equation  where  values  of  W  appearing 
as  the  unknown  on  the  left  side  also  appear  in  the  integrand.  Equa¬ 
tion  (C.2)  is  an  integral  f 'multi  where  values  of  W  appearing  in  the 
integrand  are  already  known  f rom  the  solution  of  (C.l).  The  value  of 
W  on  the  left  is  not  generally  from  the  domain  of  the  integral. 

Second,  the  function  ( 1 /r^ ) ( dv ?/Dn)  in  both  equations  is  singu¬ 
lar  at  the  point  where  r ^  -  0,  which  can  occur  only  in  the  case  of 
zero  elevation.  That  term  contributes  the  amount  required  to  maintain 


L 
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continuitv.  Sinco  the  integrals  on  a  general  surface  are  difficult, 
we  illustrate  by  considering  a  plane  earth. 

On  a  uniform  plane  surface,  dr^/fhi  =  0.  Then 


W(r  )  -  1+- 

U  l 


f  U  -r<- .  w 
l"J  rI  r2 


ik(r  +r  -r  ) 

W(Q)  c  U  <S  . 


(C.  3) 


We  wish  to  compare  this  with 


lim  2W(rn,  h)  =  lim  1 
h+0  h-*0 


J  rlri 


W(Q)  e 


ik(rl+r2_r0) 


Mi++)&] 


(C.4) 


For  receivers  above  a  plane  earth,  9r^/3n  =  -h/^*  Consistency  for 
Kqs.  (C. 3)  and  (C.4)  would  then  require 


Um  W(r  ,  h)  =  -2-  rQ 
h  v0 


M 


W(Q)  eik(rl+r2'r0) 
rlr2 


l1  +  kr2) 


2  2  2 

where  da  =  r,*?  dr,’.,  dO^,  r2  =  r2  +  h  »  and,  in  coordinates  centered  at 
tlu*  observation  point  P,  the  integration  point  Q  is  at  ^2^*  T^ie 
following  figure  shows  the  relationship  between  r^,  r^,  r^,  and  0^1 
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We  thus  evaluate 


lim  W(rn,  h) 
h-*0  U 


..  -ik  f  f  ,  .  ,  |fl  W(r2*  °2)  ik(ri+r2-ro>  h 

2tT  r0  /  /  r2  dr2  d02  "77 - °  7 

J  J  ]  r 


+  lim 
h-K) 


2tt 


W(r,  e  ) 

r '  dr;  deo  — - — ~  . 

2  2  2  r  i 

0  "0  1  r2 


// 
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We  begin  with  the  first  term  on  the  right  side  of  Eq .  (C.5) 
Note  that 


r^h  r^h 


is  a  function  that  in  the  limit  h  -+  0  tends  to  zero  for  all  r’  /  0. 

2  2  ^ 

For  r^  +  0,  the  term  (r^lO/Cr^  +  h  )  remains  finite,  so  the  entire 

integral  vanishes.  The  same  behavior  is  shown  from  a  different  per- 

2  2 

spective  by  regrouping  the  terms:  in  the  limit  h  -*  0,  h/(r ^  +  h  ) 

is  a  sharply  peaked  function  around  r^  -  0.  Thus  only  values  of  the 
integrand  for  r^  =  0  will  contribute.  But,  being  proportional  to  r^, 
the  integral  vanishes. 

The  remaining  term  contains  the  function 


r^h 


h2  +  ”2) 


2\ 3  /  2 


This  function  is  also  strong]  v  peaked  at  rj  =  0  in  the  limit  h  0 . 
However,  it  becomes  large  at  r^  =  0,  so  the  integral  does  not  vanish. 

The  rest  of  the  integrand  will  be  constant  in  the  infinitesimal  inter¬ 
val  around  =  0.  In  this  case,  we  have  and  ^(r^,  ©2)  v  W(Tq). 


Thus 


-62- 


11m  W(r  ,  h) 
IrH)  U 


lim 

h>0 


w(r0)  f'  jf 

J  J  77' 

*'o  ^0  \  2 


r'b  dr' 


2  +  h2\3/2 


But 


lim 

h-K) 


r^h  dr^ 


lim - 

h+0  (r2 


-h 


h2)i/2 


oo 


0 


lim 

h-K) 


-h 


(»2  +  0)1/2 


+ 


(0  +  h2)1/2 


=  1 


Thus 


lim  W(rn,  h)  =  W(rn) 
h~K) 


and  the  elevated  receiver  formula  tends  to  the  same  value  of  h  -> 
is  given  by  the  equation  for  the  function  on  the  surface. 
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SYMBOI^S 


A 

-> 

B 

C 

E 

F(t) 

H 

k 


-> 


n 


g 

t 


v 


W 

2  (x) 
6 
Ax 


e 


0 


K 


Vi 

o 

4> 

<P 

% 

* 


vector  tangent  to  surface  of  hill 

vector  tangent  to  surface  of  hill  independent  of  K 

-  > 

vector  normal  to  A  and  B 
electric  field 

LORAN  signal  as  function  of  time 
magnetic  field 
wave  number 

wave  number  of  center  component  of  pulse 
unit  normal  to  ground  surface 
index  of  refraction  of  ground 
time 

time  of  arrival  for  front  of  LORAN  pulse 
time  for  LORAN  pulse  to  reach  peak 
time  of  arrival  for  third-cycle  zero 
velocity  of  Gaussian  pulse 
attenuation  function  for  ground  wave 
ground  elevation  function 
skin  depth 

width  of  Gaussian  pulse 
permittivity  of  free  space 
relative  permittivity  of  conductor 
permeability 

permeability  of  free  space 

conductivity 

phase  of  W 

ground  wave  field 

ground  wave  field  if  no  ground  present 
phase  of  0 


w  =  frequency 

Wq  =  center  frequency  of  pulse 
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